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ABSTRACT 


Analysis  of  experimental  pressure  distributions  and  force  and 
moment  coefficients  for  a  9-deg  half-angle  spherically  blunted  cone  at 
M  -  is  and  «  =  0  to  20  deg  is  presented  based  on  an  ideal  gas  (y  =  1.4) 
three-dimensional  characteristics  solution  developed  by  General  Applied 
Sciences  Laboratory  (GASL).  Comparisons  are  also  made  with  pre¬ 
dictions  based  on  modified  Newtonian  theory.  Pressure  distributions 
along  the  windward  and  leeward  streamlines  were  well  predicted  by  the 
results  from  the  GASL  program  at  a  -  10  deg.  For  a  <  LOdeg,  comparison 
with  the  experimental  data  indicated  differences  attributed  to  viscous 
effects.  Fora  >  10  deg  the  numerical  results  were  affected  by  necessary 
program  modifications  on  the  leeward  side.  Radial  pressure  distribu¬ 
tions  along  the  body  surface  from  the  GASL  results  showed  the  develop¬ 
ment  of  low  pressure  "wells"  in  the  leeward  region  at  «  >  10  deg  .  Similar 
distributions  along  a  constant  entropy  surface  within  the  inviscid  shock 
layer  indicated  substantial  delay  in  the  formation  of  the  low  pressure 
"wells'’.  Cn  and  C„,  were  found  to  be  in  good  agreement  with  the  experi¬ 
mental  data  when  based  on  the  GASL  surface  pressure  distribution. 
Predicted  values  of  C*  were  about  40  percent  below  the  experimental 
data  over  the  entire  range,  a  =  o  to  20  deg,  in  agreement  with  a  previous 
Newtonian  based  analysis.  The  difference  is  attributed  to  the  viscous- 
induced  drag  increment.  Computed  slopes,  C $ n  and  ,  were  found  to 
be  well  predicted  by  the  integrated  GASL  results,  whereas  the  modified 
Newtonian  prediction  was  in  poorer1  agreement  with  the  trends  of  the 
experimental  data. 
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A  Surface  area/Rn2 
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Ca  Axial -force  coefficient,  axial  force /q^Ab 

CD  Drag  coefficient,  D/q^Ab 


^Unless  otherwise  noted,  all  lengths  are  nondimensionalized  with 
respect  to  the  nose  (sphere)  radius. 
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D 
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i 
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Rr 

Rb 


He  =o ,  i 

Rd 


Nose  drag  coefficient  (0.964  for  a  sphere-cone) 

General  force  coefficient  representing  the  total  force  acting 
on  vehicle,  total  force/q^At 

Lift  coefficient,  Lq^As 

Resultant-moment  coefficient  vector  obtained  by  summing 
moments  of  all  forces  about  point  uO".  moment  /qM  A[,  I)h Rn 

Pitching-moment  coefficient,  magnitude  of  Y-component  of 
resultant -moment  coefficient  vector 

Pitching -moment  coefficient  slope,  PC tn  ’da,  deg" 1 

Normal-force  coefficient,  normal  force  qMAb 

Normal-force  coefficient  slope,  <9C,\''c)a,  deg-1 

Pressure  coefficient,  ip  - 

Chapman-Rubesin  viscosity  coefficient,  tp..v  ;  Pc 1  Pr-'H  w~> 
Drag  force  along  the  free-stream  velocity  vector 
Vehicle  base  diameter 
Vehicle  nose  (sphere)  diameter 
Stagnation  enthalpy,  ft2 /sec 2 

Lift  force  normal  to  free-stream  velocity  vector 
Body  length,  in. 

Mach  number 

Force  normal  to  body  axis  of  symmetry 
Unit  inward  normal  vector  to  surface 
Pressure 

Free-stream  stagnation  pressure,  lbf/in.2 
Normal-shock  pitot  pressure,  lbf/in,2 
Free-stream  dynamic  pressure,  '2,  lbf/in.2 

Body  radius  in  body-fixed  axis  system 
Vehicle  base  radius 

Free-stream  unit  Reynolds  number,  in.-* 

Free-stream  Reynolds  number  based  on  body  length 

Vehicle  nose  (sphere)  radius,  ft 

Radius  vector  from  point  "O"  to  body  point 
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Entropy,  dimensionless 

Surface  distance  from  the  geometric  nose  of  the  model 
Temperature,  °K 
Free-stream  velocity,  ft/sec 

Viscous  parameter,  \]M  \/cT  /v'Re^ 

Coordinate  axis  normal  to  Z  positive  downward 

Right-hand  coordinate  axis  normal  to  x-z  plane 

Coordinate  axis  along  the  free-stream  velocity  vector  (wind- 
fixed  axis)  or  along  the  axis  of  symmetry  (body-fixed  axis) 

Cheng's  axial  distance  parameter,  [(Zb  +  1 ) /2  ]  ac2  /  v  *Ci)n 
Center  of  pressure  location 
Angle  of  attack,  deg 

Ratio  of  specific  heats  (1,40  used  herein) 

[y  -  !)/(>•  +  1) 

Angle  defined  by  cos  if  =  V  •  n  /  |  V  | 

Polar  angle  measured  in  plane  perpendicular  to  axis  of  sym¬ 
metry  from  most  windward  streamline 

Cone  half- angle,  deg 

SUBSCRIPTS 

1,2.3  Denoted  values  of  parameters  at  limits  of  integration  for  each 
region,  at  special  values  of  Z  and  0 ,  or  values  of  force  coef¬ 
ficients  for  each  like  numbered  region 

B  Body  axes 

b  Base 

i  Initial  value  in  streamline  program 

n  Nose,  or  normal  direction 

o  Stagnation  conditions 

w  Wind  axes,  or  wall 

\i ,  y,i  i  =  I3,w  ,  that  component  of  force  in  region  acting  in  \i  or  /.i 

direction 

~  Undisturbed  free-stream  conditions 
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SECTION  I 
INTRODUCTION 


One  function  of  experimental  research  is  to  test  the  applicability 
and  range  of  validity  of  theoretical  models.  In  fluid  mechanics,  as  the 
velocity  increases  many  problems  appear  such  as  compressibility,  inter¬ 
actions  of  inviscid  and  viscous  flow  fields,  and  so-called  '  real  gas  1 
effects.  On  the  one  hand,  the  problem  is  often  complicated  by  the  in¬ 
adequacy  of  the  experimental  results  to  clearly  indicate  the  relative 
magnitudes  of  separate  but  coupled  effects.  On  the  other  hand,  theo¬ 
retical  models  usually  are  only  applicable  to  simple  geometric  and  fluid 
mechanical  models.  Even  with  the  aid  of  large  digital  computers,  'exact'1 
solutions  are  not  possible  for  many  problems  of  current  interest  in  gas 
dynamics. 

In  recent  years  considerable  attention  has  been  given  to  the  com¬ 
parisons  of  ideal  gas  and  perfect  dissociating  and  ionizing  gas  character¬ 
istics  solutions  for  sphere -cones  at  zero  angle  of  attack,  but  similar 
comparisons  have  not  been  possible  at  angles  of  attack  because  of  the 
lack  of  adequate  theoretical  models  for  analysis  of  the  three-dimensional 
flow  fields.  The  three-dimensional  characteristics  theory  used  in  the 
present  study  was  developed  at  the  General  Applied  Sciences  Laboratory 
(GASL)  by  Moretti  et  al.  (Ref.  1).  The  only  known  published  results 
from  this  theory  for  axisymmetric  bodies  are  calculations  done  by 
Jenkins  (Ref.  2)  for  spherically  blunted  cones  at  0-  and  5 -deg  angles  of 
attack  over  the  range  of  altitudes  from  100,  000  to  260,  000  ft,  Mach  num¬ 
bers  from  9  to  30,  and  cone  half-angles  of  10  and  12  deg. 

The  present  report  presents  comparisons  between  an  ideal  gas, 
inviscid,  three-dimensional  characteristic  solution  and  experimental 
data  for  surface  pressure  distributions,  forces  and  moments  of  a 
0.  3  nose-to-base  radius  ratio,  9 -deg  half-angle  spherically  blunted 
cone  at  =  18  over  a  range  of  angles  of  attack  from  0  to  20  deg.  The 
experimental  data  were  taken  in  a  Reynolds  number  range  where  strong 
viscous  effects  on  zero-lift  drag  were  previously  found  by  Whitfield  and 
Griffith  (Ref.  3).  These  experimental  data  at  angle  of  attack  will,  how¬ 
ever,  be  used  to  assess  the  limitations  and  validity  of  the  characteristics 
solution. 

Since  the  wind  tunnel  experimental  data  were  taken  in  nitrogen,  it 
was  necessary  to  base  the  calculations  on  ideal  gas  (y  =  1.4)  properties 
rather  than  the  equilibrium  air  data  available  in  the  machine  program. 
Previously  published  zero-lift  experimental  data  taken  in  the  AEDC, 
von  Karman  Gas  Dynamics  Facility  (VKF),  hotshot  tunnels  by  Whitfield 
and  Griffith  (Ref.  3),  Griffith  and  Lewis  (Ref.  4),  and  Lewis  and  Whitfield 
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(Ref,  5)  were  found  to  be  in  reasonably  good  agreement  with  their  ideal 
gas  analyses. 

Most  of  the  calculations  to  be  presented  herein  were  made  using  the 
GASL  wind-fixed  axis  program  (Ref.  1).  However,  during  the  latter 
stages  of  the  calculations  a  new  version  of  the  program  was  published  by 
GASL  (Ref.  6)  and  was  used  for  some  of  the  calculations  presented.  The 
new  program,  denoted  herein  as  the  body -fixed  axis  program,  is  a  modi¬ 
fication  of  the  wind-fixed  axis  program  for  axisymmetric  bodies  such 
that  the  calculation  is  done  in  body-fixed  rather  than  wind-fixed  axes. 

This  procedure  simplifies  the  analysis  and  somewhat  improves  the  ac¬ 
curacy  of  the  calculation.  However,  a  comparison  will  be  shown  at  one 
angle  of  attack  which  indicates  that  the  differences  are  small,  and  thus 
the  large  body  of  calculations  at  all  angles  of  attack  was  not  recomputed 
using  the  new  program.  Included  in  the  new  program  is  a  procedure  for 
tracing  the  streamline  patterns  on  the  surface  and  within  the  inviscid 
shock  layer.  This  program  is  important  for  many  applications  and  was 
used  in  the  present  study  to  determine  streamline  patterns  on  the  body 
and  on  one  constant  entropy  surface  within  the  shock  layer. 

Calculations  were  made  at  a  -  5,  10,  and  20  deg  for  comparisons 
with  experimental  data  from  the  AEDC-VKF  50-in.  and  100-in.  hyper¬ 
velocity  tunnels  (Gas  Dynamic  Wind  Tunnels,  Hypersonic  (H)  and  (F)) 
and  the  48-in.  shock  tunnel  of  the  Cornell  Aei’onautical  Laboratory  (CAL), 
A  machine  program  was  developed  by  the  first  author  and  M.  R.  Reed 
which  transformed  the  wind-fixed  axis  solution  for  the  body  surface  pres¬ 
sure  distribution  to  the  corresponding  body-fixed  axis  and  interpolated 
the  resulting  pressure  distribution  for  comparison  with  the  experimental 
data.  An  additional  machine  program  was  written  to  integrate  the  result¬ 
ing  pressure  distribution  and  compute  the  normal-force  coefficient,  Cn  , 
and  pitching-moment  coefficient,  Cm  ,  for  comparison  with  available 
experimental  data.  Also,  for  comparison  with  the  numerical  results 
based  on  the  GASL  program  and  experimental  data,  the  results  of  the 
modified  Newtonian  theory  are  presented.  In  addition  to  the  above  com¬ 
parisons,  normal-force  and  pitching-moment  slopes,  Cp,a  and  C^,  were 
computed  and  are  compared  with  numerical  differentiation  of  the  experi¬ 
mental  data. 


SECTION  II 
INVESTIGATION 


2.1  THEORETICAL  CONSIDERATIONS 

The  numerical  analysis  presented  herein  used  the  GASL  wind-fixed 
axis  program  to  obtain  the  body  surface  pressure  distribution  for  a  range 
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of  angles  of  attack.  The  wind-fixed  axis  solution  was  then  transformed 
to  body-fixed  axis  and  interpolated  for  integration  and  comparison  with 
the  experimental  data.  A  brief  discussion  of  the  GASL  program  is  pre¬ 
sented  below. 

The  initiaL  supersonic  data  were  obtained  from  an  inverse  blunt 
body  solution  for  a  unit  sphere  at  M*,  =  18  in  an  ideal  gas  at  one  atmos¬ 
phere  pressure.  The  resulting  blunt  body  flow  field  solution  was  inter¬ 
polated  along  a  surface  normal  to  the  wind-fixed  axis.  The  GASL 
program  uses  the  method  of  characteristics  in  three  dimensions  to  com¬ 
pute  from  these  initial  data  the  fluid  properties  in  the  inviscid  supersonic 
flow  field  about  spherically  capped  lifting  bodies.  The  technique  is  such 
that  computations  are  made  in  two  mutually  perpendicular  planes  using 
the  intersections  of  the  characteristic  surfaces  in  the  flow  field  with 
these  planes  as  characteristic  lines.  The  intersection  of  the  character¬ 
istic  lines  represents  a  solution  point  which  is  common  to  both  planes. 

For  the  present  calculations,  a  cylindrical  coordinate  system  was  used 
with  Zw  measured  from  the  center  of  the  sphere  in  the  wind  direction, 
the  radius  H,  normal  to  the  Zw  axis,  and  the  poLar  angle,  0W,  measured 
from  the  most  windward  streamline.  For  the  GASL  calculations  pre¬ 
sented,  the  shock  layer  was  defined  by  eleven  positions  between  the  body 
and  the  shock  and  10-deg  increments  in  0W .  The  calculations  were  made 
along  the  Zw  axis  such  that  the  solutions  were  obtained  for  each  0SV  between 
0  and  180  deg  in  the  Zw  plane,  and  then  a  subsequent  Zv,  plane  was  selected 
and  the  procedure  was  repeated. 

Inasmuch  as  the  solutions  were  obtained  in  Zw  planes  with  meas¬ 
ured  about  the  Zw  axis,  it  was  necessary  to  shift  the  7.w  axis  periodically 
as  the  calculation  proceeded  along  the  body.  Whenever  the  7,w  axis 
deviated  by  more  than  a  fixed  amount  from  the  axis  of  symmetry  of  the 
body,  a  "change-of-frame'1  was  made.  Figure  1  shows  a  sketch  of  the 
cone  nomenclature  and  coordinate  system.  The  criteria  for  determining 
how  often  a  change-of-frame  was  necessary  are  arbitrary,  and  the  results 
of  different  scheduled  changes  -of-frame  are  reported  herein. 

The  transformation  of  the  GASL  wind-fixed  axis  pressure  distribu¬ 
tion  was  made  using  the  appropriate  equations  for  a  selected  body  point 
from  which  the  equivalent  Zw  and  0W  were  determined.  The  GASL  results 
were  then  interpolated  by  a  four -point  Lagrangian  interpolation  formula 
to  obtain  the  pressure  at  desired  points.  The  resulting  interpolated  pres¬ 
sure  distribution  was  integrated  by  the  trapezoidal  rule  to  obtain  the 
forces  and  moments.  Also  for  comparison,  the  modified  Newtonian 
theory  was  used  to  calculate  a  pressure  distribution  and  then  integrated 
to  obtain  the  forces  and  moments. 
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A  detailed  discussion  of  the  pertinent  equations  and  techniques  used 
in  the  present  analysis  is  given  in  the  Appendix. 

The  GASL  body-fixed  axis  three-dimensional  characteristics  solu¬ 
tion  (Ref.  6)  was  essentially  that  described  above  except  that  the  calcula¬ 
tion  was  done  in  body-fixed  rather  than  -wind-fixed  axes.  The  initial 
data  were,  however,  obtained  by  rotating  the  wind-fixed  axis  sphere 
solution  to  obtain  the  initial  data  in  a  plane  perpendicular  to  the  body 
axis  of  symmetry. 


2.2  EXPERIMENTAL  STUDIES 

The  experimental  data  presented  herein  were  obtained  in  two 
AEDC-VKF  hotshot  tunnels.  Tunnels  H  and  F,  and  the  48-in.  reflected 
shock  tunnel  of  the  Cornell  Aeronautical  Laboratory  (CAL).  The  AEDC- 
VKF  and  CAL  tests  were  conducted  in  nitrogen  and  air,  respectively. 

All  data  were  taken  on  a  0.  3  nose-to-base  radius  ratio  9-deg  half-angle 
spherically  blunted  cone  over  an  angle -of- attack  range  of  0  to  20  deg. 

A  comparison  of  some  of  the  experimental  data  presented  herein  from 
the  VKF  and  CAL  facilities  was  previously  presented  by  Edenfield  in 
Ref.  7,  and  the  reader  is  referred  to  that  report  for  a  detailed  descrip¬ 
tion  of  the  experimental  programs.  However,  tables  are  presented 
herein  which  list  all  the  experimental  data  used  in  the  present  study. 

Table  I  gives  the  pressure  distribution  data  from  the  VKF  tunnels, 
referenced  to  the  measured  test  section  pitot  pressure.  Each  column 
gives  a  gage  location  on  the  model,  and  the  coordinates  of  each  orifice 
are  noted.  The  axial  distance  is  referenced  to  the  nose  radius  and 
measured  from  the  center  of  the  spherical  cap,  whereas  the  surface 
distance  is  measured  from  the  nose  of  the  model.  The  model  radial 
angle,  *  was  measured  from  the  most  windward  streamline  direction. 
The  test  conditions  associated  with  the  data  presented  in  Table  I  are 
given  in  Table  II. 

The  data  for  the  CAL  investigation  were  obtained  from  Ref.  8.  The 
pressure  data  are  presented  in  Table  III  in  a  modified  form,  and  the 
tunnel  test  conditions  are  given  in  Table  IV.  The  data  given  in  Ref.  8 
were  referenced  to  a  computed  p'  based  on  previous  tunnel  calibrations. 
These  data  were  modified  so  that  at  zero  angle  of  attack  the  measured 
stagnation  point  pressure  was  used  as  the  reference,  and  at  angles  of 
attack  the  pressure  measured  nearest  the  stagnation  point  was  corrected 
to  the  stagnation  point  by  modified  Newtonian  theory.  The  latter  value  is 
denoted  herein  as  p^*. 

Finally,  tabulations  of  all  the  force  and  moment  coefficients  are 
presented  in  Tables  V  and  VI.  The  experimentally  determined  coefficients 
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and  the  associated  tunnel  test  conditions  are  given  in  Table  V,  and  the 
computed  coefficients  are  given  from  the  GASL  and  modified  Newtonian 
numerical  solutions  in  Table  VI. 


SECTION  III 

RESULTS  AND  DISCUSSION 


3.1  ZERO-LIFT  PRESSURE  DISTRIBUTIONS 

An  obvious  test  of  any  three-dimensional  flow  field  solution  is  com¬ 
parison  with  a  two-dimensional  or  axisymmetric  solution  where  both 
theories  are  applicable,  viz.,  at  zero  angle  of  attack.  Such  a  com¬ 
parison  is  shown  in  Fig.  2  between  the  GASL  solution  and  the  results 
from  an  AEDC  axisymmetric  ideal  gas  characteristics  program  devel¬ 
oped  by  ARO,  Inc.  The  comparison  is  made  in  terms  of  Cp/2 ac*  against 
Zc  =  |  (Zn  0/2]  Oc1'  v'  fC[)n  .  The  correlation  parameter,  Zc  ,  was  pro¬ 
posed  by  Cheng  (Ref,  9),  and  Griffith  and  Lewis  (Ref.  4)  modified 
Cheng's  correlation  parameter,  p/MJ' y °c  ,  to  read  Cp/2ctc2  in  order 
to  account  for  the  fact  that  the  condition  imposed  by  Cheng,  p  >  >  p^  , 
frequently  is  not  satisfied.  An  improved  correlation  with  experimental 
data  was  obtained  in  this  manner. 

Some  experimental  data  from  AEDC-VKF  and  CAL  are  also  shown 
in  Fig,  2.  Of  primary  interest  here  is  the  comparison  between  the 
AEDC  ideal  gas  axisymmetric  solution  and  the  GASL  wind- fixed  axis 
solution.  The  differences  are  believed  to  be  in  the  computation  tech¬ 
niques  used  and  the  numerical  accuracy  of  each  solution.  As  noted 
earlier,  the  GASL  program  obtains  a  solution  at  each  point  represented 
by  the  intersection  of  two  characteristic  lines  in  a  <9W  plane.  This  point 
also  lies  in  a  Zw  plane,  and  the  points  between  the  body  and  the  shock  in 
each  plane  are  interpolated  at  eleven  points  in  the  R  direction  and 
eighteen  points  in  the  9  direction.  At  angles  of  attack,  large  variations 
in  the  slopes  of  the  characteristics  lines  occur  from  the  windward  to 
the  leeward  side  of  the  body,  and  this  implies  large  deviations  in  their 
intersections  from  a  single  Zw  plane.  This  requires  the  convergence 
criteria  to  be  relatively  loose  at  angle  of  attack,  and  thus  the  accuracy 
of  the  numerical  solution  suffers  at  zero  angle  of  attack.  However, 
when  the  convergence  criteria  in  the  AEDC  axisymmetric  program  are 
not  met,  the  mesh  size  is  reduced  by  doubling  the  number  of  points  along 
a  right  running  characteristic  from  the  shock  to  the  body.  This  latter 
procedure  generates  a  finer  net-work  of  points,  and  the  solution  is  in¬ 
herently  more  accurate.  Therefore,  the  difference  betwreen  the  AEDC 
and  GASL  wind-fixed  axis  solution  at  zero  angle  of  attack  must  be  viewed 
as  an  indication  of  the  built-in  inaccuracy  within  the  GASL  method. 
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The  differences  between  the  experimental  data  and  the  AEDC  axisym 
metric  characteristics  solution  shown  in  Fig.  2  have  recently  been  in¬ 
vestigated  by  Lewis  and  Whitfield  (Ref.  5).  The  results  of  that  study 
showed  that  the  difference  can  be  accounted  for  on  the  basis  of  a  com¬ 
bined  inviscid- viscous  interaction  solution.  They  showed  that  for  the 
range  of  conditions  of  interest  in  the  present  study,  the  correlation  of 
Cp/2 oc2  against  Zc  as  given  previously  by  Griffith  and  Lewis  (Ref.  4)  was 
in  good  agreement  Avith  the  results  of  the  inviscid- viscous  interaction 
solution.  Therefore,  at  zero  angle  of  attack  the  differences  between  the 
experimental  data  and  the  AEDC  inviscid  characteristics  solution  are 
believed  to  be  primarily  caused  by  viscous  effects,  and  the  differences 
between  the  two  inviscid  solutions  are  believed  to  be  numerical. 

When  considering  the  angle -of-attack  results,  one  should  keep  in 
mind  these  comments  which  only  strictly  apply  at  zero  angle  of  attack. 
The  effects  of  numerical  accuracy  and  \dscous  interaction  will  remain 
important  at  angles  of  attack  and  particularly  at  small  angles  of  attack. 
Effects  of  viscous  interaction,  at  least  on  the  windward  side,  will  be 
reduced  with  increasing  angle  of  attack. 


3.2  SURFACE  PRESSURE  DISTRIBUTIONS  AT  ANGLES  OF  ATTACK 


3.2.1  Ehects  of  Change-of-Frame 


In  performing  the  calculations  with  the  GASL  wind-fixed  axis  pro¬ 
gram  at  angles  of  attack,  it  was  necessary  to  select  positions  along  the 
body  at  which  changes -of-frame  of  reference  were  made.  Such  changes - 
of-frame  were  made  by  specifying  the  Zw  at  ■which  the  change  was  to 
occur  and  the  incremental  change  in  Xw .  Since  both  quantities  Avere 
arbitrarily  chosen,  it  was  of  interest  to  see  the  effects  of  these  quanti¬ 
ties  on  the  solution.  To  investigate  these  effects  several  computer 
solutions  were  obtained  for  the  10-deg  angle -of-attack  case  with  a  wide 
variation  in  the  change-of-frame  schedules.  The  results  of  these  com¬ 
parisons  will  now  be  presented. 

Three  different  solutions  were  obtained  at  an  angle  of  attack  of 
10  deg  with  three  change-of-frame  schedules.  Figure  3  is  a  sketch  of 
the  model  showing  the  changes -of-frame  for  the  three  solutions  denoted 
solutions  1,  2,  and  3.  Solutions  1  and  2  differ  only  after  the  fourth 
change-of-frame,  and  solution  number  3  was  selected  such  that  AXW  =  0.2 
that  is,  whenever  A  Zw  tan  a  =  0.2  a  change-of-frame  was  made.  All  three 
solutions  were  found  to  be  practically  identical  to  change-of-frame  num¬ 
ber  5  for  solution  number  2  and  to  change-of-frame  number  10  for  solu¬ 
tion  number  3.  However,  solution  number  1  departed  radically  from 
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solutions  number  2  and  3  further  downstream.  Figure  3b  shows  the 
radial  surface  pressure  distribution  for  the  three  solutions  at  two  Zb 
stations  where  experimental  data  were  available  for  comparison.  The 
comparison  shows  solution  number  2  is  considerably  improved  over 
solution  number  1.  It  is  interesting  to  note,  however,  that  a  further 
increase  in  the  number  of  changes -of- frame  had  no  significant  effect 
on  the  distribution  as  shown  by  solution  number  3.  Since  the  criteria 
for  making  a  change -of -frame  are  arbitrary,  it  does  appear  that  one 
should  investigate  the  effects  by  making  comparisons  of  the  type  shown 
above. 

The  solution  from  the  GASL  body-fixed  axis  program  is  also  shown 
in  Fig.  3b.  Comparison  of  the  body -fixed  axis  solution  and  solution 
number  3  shows  excellent  agreement  except  near  the  most  windward 
and  most  leeward  streamlines.  There  are  also  some  small  differences 
in  the  region  of  minimum  pressure.  It  is  also  of  interest  to  note  that 
the  experimental  pressure  distributions  along  the  most  windward  (#b  =  0) 
and  most  leeward  =  ~)  streamlines  are  in  better  agreement  with  the 
body- fixed  axis  solution.  Because  of  this  and  certain  overall  improve¬ 
ments  in  smoothness  of  the  solution,  the  body -fixed  axis  solution  is 
preferred.  However,  as  noted  above,  most  of  the  calculations  for  this 
study  were  completed  before  the  body-fixed  axis  program  became  avail¬ 
able. 

Solutions  were  obtained  for  a  -  3,  10,  and  20  deg  with  change-of-frame 
schedules  similar  to  those  of  solution  number  1  shown  in  Fig.  3a.  Only 
after  examination  of  solution  number  1  and  the  ’’wells  1  of  low  pressure 
noted  were  the  two  other  solutions  obtained  for  a  =  10  deg.  The  «  =  5-deg 
case  did  not  exhibit  the  formation  of  any  low  pressure  wells.  However, 
for  a  =  20  deg  the  problem  of  low  pressure  wells  became  more  acute,  and 
the  solution  failed  to  converge  on  the  leeward  side  at  a  Z-A  station  only 
about  one-third  of  the  total  length  of  the  body.  In  order  to  continue  the 
solution  over  the  remainder  of  the  body,  it  was  necessary  to  delete  0W 
planes  in  the  leeward  region  downstream  of  the  Zw  plane  where  the  failure 
occurred.  This  was  accomplished  by  limiting  the  maximum  value  of  0W 
to  less  than  the  value  where  the  failure  occurred.  Before  the  a  -  20-deg 
solution  was  continued  an  a  =  15-deg  solution  was  obtained  with  a  change- 
of-frame  criteria  similar  to  that  in  solution  3  above,  viz.  AZW  tan  a  =  0.2. 
Difficulties  similar  to  those  encountered  in  the  a  -  20-deg  solution  were 
observed  for  the  «  =  15-deg  case.  Except  in  the  latter  case,  the  failure 
occurred  about  two-thirds  of  Znmax-  However,  it  was  necessary  again 
to  delete  6W  planes  to  continue  the  solution.  Because  of  this  result,  it 
was  felt  that  increasing  the  changes -of-frame  in  the  20-deg  case  would 
not  improve  the  solution  significantly,  since  deleting  planes  in  the  lee¬ 
ward  region  is  tantamount  to  removing  the  leeward  boundary  condition 
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from  the  solution,  viz.  ,  the  requirement  that  the  cross  velocity  be  zero 
along  the  most  leeward  streamline.  Little  can  be  said  regarding  the 
validity  of  solutions  obtained  under  conditions  where  6n  planes  are 
deleted.  Therefore,  only  representative  radial  distributions  will  be 
presented  to  illustrate  the  character  of  the  solutions  upstream  and  down¬ 
stream  of  the  Zw  plane  where  the  solutions  failed  to  converge.  Figure  4 
shows  the  change -of -frame  schedules  for  all  angles  of  attack  used  in 
the  wind -fixed  axis  solutions. 

3.2.2  Windward  and  Leeward  Streamline  Pressure  Distributions 


The  computed  longitudinal  pressure  distributions  on  the  windward 
(#3  =  o)  and  leeward  =  rr)  streamlines  on  the  conical  portion  of  the 
body  are  compared  with  experimental  data  in  Figs.  5a,  b,  and  c  for 
a  =  5,  10,  and  20  deg.  Lack  of  experimental  data  at  a  =  15  deg  prevented 
a  similar  comparison  for  that  case.  The  experimental  data  show  con¬ 
siderable  scatter,  but  in  comparison  with  the  computed  distributions, 
the  agreement  is  reasonably  good.  The  agreement  is  improved  with  in¬ 
creasing  angle  of  attack,  especially  along  the  most  windward  streamline, 
and  as  noted  above,  viscous  effects  are  reduced  with  increasing  angle 
of  attack.  At  a  -  10  deg  the  distributions  along  the  most  windward  and 
leeward  streamlines  from  the  body-fixed  solution  are  also  shown,  and 
the  results  are  seen  to  be  in  better  agreement  than  the  corresponding 
wind-fixed  axis  solutions. 

At  a  =  5  and  10  deg  along  the  leeward  streamlines  the  agreement 
between  the  experimental  and  numerical  pressure  distributions  is  not  as 
good  as  for  the  windward  streamlines.  The  agreement  improves  some¬ 
what  with  angle  of  attack.  However,  because  of  the  low  pressure  region, 
experimentally  there  will  be  a  thick  laminar  boundary  layer  interacting 
with  the  inviscid  outer  flow.  This  viscous  interaction  will  induce  a  pres¬ 
sure  increment  as  noted  previously  for  zero  angle  of  attack.  Since  the 
solution  along  the  leeward  streamline  at  20-deg  angle  of  attack  only 
existed  over  approximately  the  first  one-third  of  the  body,  the  com¬ 
parison  shown  in  Fig.  5c  is  interesting,  but  conclusive  statements 
regarding  the  results  are  not  possible. 

It  is  well  known  that  the  nose-dominated  region  of  spherically 
blunted  slender  cones  extends  many  nose  diameters  downstream.  That 
is,  for  a  cone  of  small  half-angle  {e.  g.  ac  <  9  deg  ),  the  inviscid  blunt 
cone  surface  pressure  distribution  will  require  more  than  60  nose  radii 
before  approaching  the  inviscid  sharp  cone  value.  As  the  cone  angle  in¬ 
creases  the  nose -dominated  region  is  reduced  in  extent.  The  pressure 
distribution  is  characterized  by  a  strong  overexpans  ion  and  recompres¬ 
sion  to  the  sharp  cone  value.  It  is  of  interest  to  note  at  zero-lift  that 
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the  cone  under  investigation  ends  at  approximately  the  minimum  in  the 
overexpansion  and  recompression  region.  It  is  also  of  interest  to  note 
from  Figs.  5a  and  b  that  the  same  effects  of  the  nose -dominated  region 
are  observed  along  the  windward  streamlines  at  a  =  5  and  10  deg. 

Axisymmetric  characteristic  solutions  were  obtained  for  spherically 
blunted  cones  with  half-angles  crc  =  9  +  a  deg  .  Figure  5d  shows  a  com¬ 
parison  of  the  axisymmetric  solutions  for  =  14  and  19  deg  and  the 
11  corresponding"  GASL  solutions  for  a  =  5  and  10  deg.  Of  primary  interest 
here  is  the  approximately  constant  difference  between  the  solutions  for 
the  range  of  Zw  shown  and  the  smaller  difference  between  solutions  as  ac 
or  a  increases.  This  is  in  agreement  with  the  results  and  comments 
given  for  a  =  o.  It  is  also  of  interest  to  note  the  small  effects  of  cross 
flow  along  the  most  windward  streamline  as  indicated  by  the  good  agree¬ 
ment  between  the  axisymmetric  and  three-dimensional  solutions.  It 
thus  appears  that  an  axisymmetric  characteristics  solution  adequately 
predicts  the  most  windward  streamline  pressure  level  and  the  over¬ 
expansion  and  recompression  regions  for  an  equivalent  angle,  oCe,  y  = 
ctc  +  a.  Except  in  the  asymptotic  limit  of  large  Zb  or  large  cone  angle, 
oc ,  modified  Newtonian  is  unable  to  predict  the  correct  pressure  level 
and,  in  any  event,  is  unable  to  predict  the  overexpansion  and  recom¬ 
pression  pressure  region.  The  axisymmetric  calculations  gave 
p/p^  =  0.0594,  0.1230,  and  0.2440  fora  =  5,10,  and  20  deg,  respectively,  for 
the  limiting  surface  pressure  ratio,  and  modified  Newdonian  theory  gave 
corresponding  values  of  0.057,  0.  105,  and  0.  236.  A  comparison  with 
the  experimental  data  of  Table  I  shows  good  agreement  with  these  values. 
The  limits  of  applicability  of  the  Newtonian  theory  will  be  clearer  from 
comparisons  of  predicted  and  experimental  radial  pressure  distributions 
at  angles  of  attack. 

3.2.3  Radial  Pressure  Distributions 


Comparisons  of  the  radial  pressure  distributions  at  several  posi¬ 
tions  along  the  body  computed  by  the  GASL  wind-fixed  axis  program  and 
modified  Newtonian  are  presented  in  Fig.  6  for  a  =  5  deg .  The  available 
experimental  data  are  also  shown  for  comparison.  Only  at  =  10.675 
were  data  available  for  a  range  of  the  viscous  hypersonic  parameter,  . 
The  trends  of  the  distributions  from  the  GASL  wind -fixed  axis  solution 
and  experiment  are  in  good  agreement.  As  noted  before,  the  experi¬ 
mental  data  are  strongly  viscous  affected  especially  at  this  low  angle  of 
attack.  It  should  be  noted  here  that  the  CAL  experimental  data  are  for 
M  -  13  ,  whereas  the  calculation  was  made  for  =  18  .  It  would  be 
surprising  indeed  if  the  differences  shown  were  strongly  Mach  number 
dependent.  Moreover,  the  Reynolds  numbers  for  the  experiments  were 
sufficiently  high  (v^  =  0.03  v)  that  it  would  be  equally  surprising  if  the 
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differences  between  experiment  and  numerical  results  were  attributable 
entirely  to  viscous  interaction.  Thus,  further  calculations  and  experi¬ 
ments  at  the  CAL  conditions  would  be  of  interest  at  a  =  5  deg  . 

Modified  Newtonian  theory  is  also  shown  in  Fig.  6  for  comparison 
with  the  GASL  solution  and  experimental  data.  The  comparison  shows 
reasonably  good  agreement  over  the  windward  side  and  far  downstream 
of  the  nose.  It  is  under  these  conditions,  of  course,  where  the  shock 
wave  is  nearly  parallel  to  the  body  that  one  would  expect  Newtonian 
theory  to  better  apply;  however,  even  over  the  windward  side  the  GASL 
solution  better  predicts  the  trends  of  the  experimental  data. 

Radial  distributions  for  a  =  10  deg  are  shown  in  Fig.  7  for  compari¬ 
son  with  the  considerably  larger  body  of  experimental  data.  It  was,  in 
fact,  this  body  of  data  which  influenced  the  decision  to  obtain  the  GASL 
body-fixed  axis  solution  for  this  condition.  In  contrast  to  the  a  -  5-deg 
experimental  data,  the  experimental  scatter  here  appears  to  be  con¬ 
siderably  reduced,  and  the  agreement  with  the  numerical  surface  pres¬ 
sure  distribution  is  quite  good  to  Or  =  120  deg.  The  largest  body  of 
experimental  data  exists  for  this  case  at  Zb  =  10.675,  and  as  in  the 
a  -  5-deg  case,  there  was  no  consistent  trend  of  the  experimental  data 
with  Reynolds  number.  Modified  Nev/tonian  theory  is  seen  to  be  in 
better  agreement  with  the  experimental  data  in  the  windward  region  than 
in  the  previous  a  ■  5  deg  case.  This  is  not  surprising  since  the  shock 
wave  lies  closer  to  the  body  along  the  windward  side  and  is  more  nearly 
parallel  with  the  body.  Modified  Newtonian  theory  also  again  fails  in  the 
leeward  region,  and  comparisons  of  the  GASL  solution  and  experimental 
data  also  show  poor  agreement  on  the  leeward  side  of  the  body. 

The  discovery  of  the  inviscid  low'  pressure  "wells"  on  the  leeward 
side  of  the  body  at  a  >  10  deg  and  their  apparent  strong  effects  on  the 
solution  was  one  of  the  more  interesting  aspects  of  the  present  study. 
These  low  pressure  regions  have  been  observed  experimentally  on  other 
axisymmetric  bodies  at  low  Mach  numbers  and  high  Reynolds  numbers. 
However,  the  experimental  data  at  Mm  -  18  did  not  indicate  their  existence. 
A  well  was  first  noted  in  the  distributions  shown  here  at  Zb  =  5.250.  How¬ 
ever,  it  was  found  that  once  the  well  began  to  develop  the  minimum  pres¬ 
sure  tended  to  drop  rather  rapidly  with  increasing  Zb  and  ultimately  led 
to  the  failure  of  the  solution  to  converge.  The  a  =  10-deg  calculation  was 
continued,  and  the  solution  failed  to  converge  at  Zb  «  22.  The  experi¬ 
mental  data  do  not  appear  to  support  the  development  of  the  low  pres¬ 
sure  wells.  However,  viscous  effects  could  mask  the  effect  and  almost 
certainly  would  strongly  affect  the  pressure  distribution  in  this  region. 

There  appear  to  be  two  simple  explanations  for  the  development  of 
the  inviscid  wells.  First,  consider  the  flow  along  a  streamline  on  the 
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body  surface.  The  fluid  expands  along  the  streamline  from  the  high 
pressure  windward  side.  However,  as  the  leeward  side  of  the  body  is 
approached  it  is  necessary  for  the  streamlines  to  be  turned  back  toward 
the  free-stream  direction  since  the  streamlines  cannot  cross.  This 
bunching  of  streamlines  near  the  most  leeward  streamline*  will  create 
a  ridge  of  high  pressure.  The  second  explanation  is  related  to  the 
alleged  loose  convergent  criterion  discussed  earlier.  At  the  higher 
angles  of  attack  the  right  running  characteristics  on  the  leeward  side 
of  the  body  are  more  nearly  parallel  with  the  body  so  that  the  distance 
between  the  intersections  becomes  large.  Therefore,  the  possibility  of 
errors  affecting  the  solution  is  increased.  The  first  explanation  appears 
physically  reasonable.  However,  in  any  numerical  solution  as  complex 
as  this  the  effects  of  numerical  inaccuracies  can  never  by  entirely 
neglected.  In  any  event,  let  it  suffice  to  say  that  a  problem  area  exists 
in  the  simple,  straightforward  application  of  the  numerical  results. 

Finally,  distributions  for  a  =  15  and  20  deg  are  shown  in  Figs.  8a  and  b, 
respectively,  to  further  illustrate  the  development  of  the  low  pressure 
wells  and  to  show  the  distributions  at  Z,\\  -=  10.675  obtained  downstream  of 
the  leeward  region  cut-off  plane.  Notice  the  excellent  agreement  between 
the  GASL  solution  and  modified  Newtonian  theory.  However,  in  view  of 
the  procedure  used  to  obtain  the  GASL  solution  (viz,,  deleting  <5W  planes) 
the  results  must  be  viewed  with  caution. 

Beyond  the  interest  in  the  comparisons  with  experimental  pressure 
distributions  shown  here,  the  results  should  be  useful  in  conjunction 
with  a  three-dimensional  boundary-layer  analysis.  Examination  of  the 
surface  pressure  distributions  for  a  =  10  deg  leads  one  to  consider  the 
effects  of  the  low  pressure  wells  and  strong  adverse  pressure  gradients 
on  a  boundary-layer  solution.  It  is  well  known  that  laminar  boundary 
layers  cannot  sustain  strong  adverse  pressure  gradients,  and  the  ad¬ 
verse  gradient  shown  for  a  =  10  deg  and  Zb  =  12.55  in  Fig.  7  would  be  suf¬ 
ficient  to  separate  a  laminar  boundary  layer.  To  further  investigate  this 
problem,  the  following  hypothetical  case  was  studied. 


3.3  PRESSURE  AND  STREAMLINE  DISTRIBUTIONS  ON  A  CONSTANT  ENTROPY  SURFACE 

To  investigate  the  behavior  of  the  pressure  distribution  within  the 
shock  layer,  a  constant  entropy  surface  (S/R  =  15.754)  was  considered  at 
a  =  10  deg  .  The  entropy  at  the  intersection  of  the  initial  surface  and  the 
bow  shock  wave  at  6>b  -  0  was  chosen.  The  initial  surface  was  a  plane 


*The  authors  are  indebted  to  E.  A.  Sanlorenzo,  General  Applied 
Sciences  Laboratory,  for  an  interesting  discussion  of  this  point. 
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normal  to  the  body  axis  at  Zg  =  -0.33787  and  intersected  the  bow  shock 
wave  and  5b  =  0  at  60.  25  deg  (acute  angle  measured  at  the  center  of  the 
sphere  from  Zw  axis).  This  value  of  the  entropy  is  the  minimum  for  the. 
complete  three-dimensional  characteristics  flow  field  solution  and  is 
independent  of  angle  of  attack  for  these  initial  surface  data.  It  should 
be  emphasized  that  this  constant  entropy  surface  (bounding  an  "entropy 
layer")  is  not  the  "edge"  of  the  boundary  layer  but  is  rather  a  surface 
within  the  inviscid  shock  Layer  along  which  certain  data  were  obtained. 

The  GASL  streamline  tracing  program  was  used  to  interpolate  the 
inviscid  shock  layer,  establish  the  coordinates  of  the  constant  entropy 
surface,  trace  the  streamlines  on  that  surface,  and  interpolate  for  the 
flow  field  variables  along  the  streamlines.  Figure  9  shows  the  cross 
sections  of  the  body,  shock  wave,  and  constant  entropy  surface  at 
selected  stations  for  a  =  10  deg.  Figure  10  shows  the  streamline  distri¬ 
butions  on  the  body  and  on  the  constant  entropy  surface.  Note  the  signif¬ 
icantly  reduced  curvature  of  the  streamlines  on  the  constant  entropy 
surface  relative  to  those  on  the  body  surface.  More  significant  com¬ 
parisons  are  shown,  however,  in  Fig.  11  where  pressure  contours  on 
the  body  surface  and  the  constant  entropy  surface  are  shown.  The 
development  of  the  low  pressure  wells  on  the  body  surface  is  easily 
seen  in  Fig.  11a,  and  a  low  pressure  well  is  seen  to  form  at  Zg  =  18 
on  the  constant  entropy  surface.  Note,  however,  that  the  wells  begin 
to  form  on  the  constant  entropy  surface  approximately  three  times 
further  downstream  than  on  the  body  surface.  It  should  also  be  men¬ 
tioned  here  that  the  solution  failed  to  converge  at  Zr  *  22  on  the  body 
surface.  Thus,  the  formation  of  the  well  and  the  associated  adverse 
pressure  gradient  on  the  constant  entropy  surface  would  indicate  that  a 
much  smaller  adverse  pressure  gradient  exists  within  the  shock  layer 
than  on  the  body  surface.  These  results,  although  only  qualitative,  are 
believed  to  be  significant  for  they  clearly  show  that  whereas  a  boundary- 
layer  solution  would  separate  based  on  the  inviscid  wall  pressure  distri¬ 
bution,  a  similar  calculation  based  on  the  boundary-layer-edge  pressure 
distribution  (simulated  here  by  the  constant  entropy  surface)  would  be 
less  subject  to  strong  adverse  pressure  gradients. 

Comparisons  between  the  experimental  surface  pressures  and  the 
computed  pressures  on  the  constant  entropy  surface  at  a  =  10  deg  are 
shown  in  Fig.  7,  The  agreement  between  experiment  and  the  constant 
entropy  surface  data  is  better  than  a  similar  comparison  with  the  pres¬ 
sure  distribution  along  the  inviscid  wall.  Because  of  the  arbitrariness 
of  the  constant  entropy  surface,  only  the  trends  should  be  considered 
here,  and  those  trends  are  in  good  agreement  with  the  experimental 
data. 
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3.4  FORCE  AND  MOMENT  COEFFICIENTS 

The  second  area  of  primary  interest  in  the  present  report  is  a  com¬ 
parison  between  calculated  and  experimentally  measured  force  and 
moment  coefficients.  Numerical  results  will  be  presented  based  on 
surface  pressure  distributions  obtained  from  the  GASL  solutions  and 
modified  Newtonian  theory.  In  addition,  numerical  and  experimental 
normal -force  and  pitching-moment  coefficients  will  be  used  to  calculate 
the  slopes,  CNa  and  Cma.  However,  before  discussing  the  results  at 
angle  of  attack  a  few  comments  will  be  made  on  previous  studies  at 
zero -lift. 

3.4.1  Zero-Lift  Drag  Studies 

The  effects  of  viscous  interaction  on  wall  measurable  quantities 
have  received  much  attention  in  VKF  in  recent  years.  The  primary 
shape  considered  in  those  studies  was  also  a  0.  3  bluntness  9-deg  half- 
angle  cone.  The  effects  of  viscous  interaction  on  zero-lift  drag  were 
studied  by  Whitfield  and  Griffith  (Refs.  3  and  10).  More  recently,  Lewis 
and  Whitfield  (Ref.  5)  re-examined  the  experimental  data  and  earlier 
analysis  at  zero-lift.  Two  important  results  from  these  studies  should 
be  noted  here:  First,  for  the  range  of  experimental  conditions  of  the 
present  study  the  total  viscous  drag  was  approximately  twice  the  inviscid 
pressure  drag.  Secondly,  the  application  of  the  most  exact  currently 
available  theories  (viz.  characteristics  and  nonsimilar  boundary-layer 
solutions)  at  zero-lift  underestimated  the  total  drag,  thus  implying  sig¬ 
nificant  contributions  by  other  so-called  "second-order"  boundary -layer 
effects  (especially  shock  generated  external  vorticity,  cf.  Ref.  5). 

3.4.2  Angle -of-Attock  Results 

Figure  12  shows  a  comparison  of  experimental  and  theoretically 
predicted  force  and  moment  coefficients  for  a  9-deg  half-angle  spher¬ 
ically  blunted  cone  over  a  range  of  angles  of  attack,  a  =  0  to  20  deg.  The 
experimental  data  shown  from  AEDC-VKF  and  CAL  were  all  chosen  at 
nominal  values  of  =  0.05  and  wall-to-stagnation  temperature  ratio 
Tw/T0  =  0.1.  The  present  results  (in  agreement  with  the  previous  results 
of  Whitfield  and  Wolny  (Ref.  11))  show  that  the  inviscid  normal-force  and 
pitching -moment  coefficients,  C\  and  Cm ,  are  in  reasonably  good  agree¬ 
ment  with  the  experimental  data  over  the  entire  angle-of-attack  range. 
However,  in  agreement  with  the  zero-lift  data,  the  axial-force  coeffi¬ 
cient,  C,\>  is  strongly  viscous  affected  over  the  angle-of-attack  range 
considered  herein.  The  trends  of  the  numerical  and  experimental 
results  are,  however,  in  reasonably  good  agreement.  For  comparison, 
Harris'  results  from  the  NASA  helium  tunnel  (Ref.  12)  and  a  R„/Rb  -  0.255 
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10-deg  half-angle  cone  at  Mm  -  19  and  Re^/in.  =  200,000  (corresponding  to 
a  =  0.016)  are  also  shown.  The  agreement  between  the  present 
analysis  for  y  =  7/5  and  the  experimental  data  in  helium  (y  =  5/3)  is 
excellent  for  a  <  7.5  deg .  *  At  a  Reynolds  number  almost  one  order  of 
magnitude  higher  than  the  AEDC-VKF  and  CAL  experimental  data,  the 
viscous -induced  drag  should,  of  course,  be  substantially  reduced.  The 
strong  effects  of  viscous  interaction  on  the  axial -force  coefficient  are 
also  shown  on  Fig.  12d  with  a  plot  of  lift-to-drag  ratio  against  angle  of 
attack.  From  this  figure  it  can  be  seen  that  the  integrated  GASL  results 
substantially  overpredict  the  lift-to-drag  ratio  because  of  the  rather 
large  viscous -induced  axial-force  increment.  Again,  it  is  also  interest¬ 
ing  to  note  that  the  NASA  helium  tunnel  data  are  in  substantial  agreement 
with  the  integrated  GASL  results. 

Modified  Newtonian  theory  was  also  used  to  compute  the  force  and 
moment  coefficients  for  the  range  of  angles  of  attack,  a  <  20 deg,  Fig.  12. 
Based  on  the  comparison  with  the  experimental  data,  the  theory  is  not 
clearly  inferior  to  the  GASL  prediction  and  appears  to  be  a  reasonably 
good  approximation  for  the  normal-force  and  pitching -moment  coeffi¬ 
cients,  Cn  and  Cm  .  Moreover,  only  relatively  small  differences  exist 
between  theoretical  predictions  of  Ca  based  on  modified  Newtonian  theory 
and  the  integrated  GASL  results.  The  area  of  largest  difference  between 
the  two  theories  occurs  at  the  lower  angles  of  attack  (a  <  10 deg)  ,  The 
significant  differences  between  the  predicted  values  of  the  normal-force 
coefficient,  Cn  ,  are  reflected  in  the  lift-to-drag  ratio  curve,  where 
again  it  should  be  noted  that  better  agreement  exists  between  the  inte¬ 
grated  GASL  results  and  the  experimental  data  than  was  found  using  the 
modified  Newtonian  theory. 

3.4.3  Stability  Derivatives 

Stability  analyses  make  use  of  normal-force  and  pitching-moment 
coefficient  slopes,  CNa  and  Cm<7.  It  was,  therefore,  of  interest  to  examine 
the  integrated  GASL  and  modified  Newtonian  results  for  predictions  of 
those  slopes.  The  coefficients  Cn  and  C,*  are  again  shown  in  Fig.  13 
based  on  the  two  theoretical  models.  In  addition,  the  variation  of  CNa 
and  Cma  with  angle  of  attack  is  also  shown.  An  experimental  uncertainty 
of  ±10  percent  was  associated  with  the  experimental  AEDC-VKF  and  CAL 
data  and  is  also  shown.  The  minimum  and  maximum  of  the  experimental 


*The  differences  in  cone  geometry,  Mach  number,  and  y  make  a 
quantitative  comparison  between  the  numerical  results  at  y  =  1.4  and  the 
helium  tunnel  experimental  data  impossible;  the  trends  of  the  experi¬ 
mental  data  are  as  expected. 
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data,  treated  as  two  separate  curves,  were  interpolated  at  2.  5-deg 
intervals  with  a  four -point  Lagrangian  interpolation  formula  over  the 
range  a  =  0  to  15  deg.  The  results  were  differentiated  numerically  with 
a  cubic  spline,  and  the  resulting  slopes  are  shown  as  shaded  areas  in 
the  figure.  Of  most  importance  here  is  the  trend  of  the  results  rather 
than  the  absolute  values  of  the  normal-force  and  pitching-moment  slopes. 
Absolute  differences  between  experiment  and  theory  are  difficult  to 
determine  because  of  the  uncertainty  in  the  experimental  data  and  the 
errors  induced  by  numerical  differentiation.  It  is  clear  from  the  results 
shown  that  the  integrated  GASL  results  over  the  range  a  =■  oto  10  deg  give 
trends  that  are  in  good  agreement  with  the  experimental  data,  whereas 
the  modified  Newtonian  theory  results  are  in  poor  agreement  with  the 
experimental  trends. 


SECTION  IV 

CONCLUDING  REMARKS 


Comparisons  have  been  made  between  an  ideal  gas  (y  =  1.4)  three- 
dimensional  characteristics  solution  (GASL),  the  modified  Newtonian 
theory  and  available  experimental  data  for  pressure  distributions,  force 
and  moment  coefficients,  and  stability  derivatives  for  a  9-deg  spher¬ 
ically  blunted  cone  over  the  angle-of-attack  range  from  0  to  20  deg.  The 
results  were  as  follows: 

1,  The  effects  of  arbitrarily  prescribed  change -of-frame  criteria 
in  the  wind-fixed  axis  solution  wrere  investigated.  The  criteria 
AZW  tan  a  =  0.2  gave  good  results  over  the  range  of  angles  of 
attack  and  axial  distance  considered  in  the  present  study.  Some 
improvement  was  found  by  using  the  body-fixed  axis  solution 
and  is  to  be  preferred  in  further  studies  of  axisymmetric  bodies 
at  angles  of  attack.  However,  the  differences  found  in  surface 
pressure  distributions  were  small,  and  the  wind-fixed  axis 
results  given  in  the  present  report  were  only  slightly  affected. 

2.  A  comparison  between  the  GASL  results  and  an  axisymmetric 
characteristic  solution  at  zero  angle  of  attack  indicated  differ¬ 
ences  in  the  pressure  distribution  along  the  body  of  about 

10  percent.  The  GASL  result  was  found  to  be  uniformly  below 
an  AEDC  ideal  gas  axisymmetric  solution.  Previous  studies 
have  shown  that  differences  between  inviscid  solutions  and 
experimental  data  of  the  order  of  10  percent  can  be  explained 
as  the  viscous -induced  pressure  increment.  Therefore,  caution 
should  be  exerted  when  using  the  GASL  wind-fixed  solution  to 
predict  pressure  distribution  at  zero  angle  of  attack. 
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3.  Predictions  of  windward  streamline  pressure  distributions 
from  the  GASL  results  were  found  to  be  in  better  agreement 
with  the  experimental  data  with  increasing  angle  of  attack,  and 
at  10-deg  angle  of  attack  the  comparison  was  quite  good.  Be¬ 
cause  of  the  necessity  of  dropping  6W  planes  at  a  =  15  and  20  deg  , 
the  comparison  of  windward  streamline  pressure  distribution 

at  a  =  20  dog  was  interesting,  but  conclusive  statements  regard¬ 
ing  the  comparison  were  not  possible. 

4.  Comparison  of  experimental  and  theoretically  predicted  radial 
pressure  distributions  based  on  the  GASL  and  modified 
Newtonian  theories  showed  the  GASL  results  to  be  uniformly 
superior  in  predicting  the  trends  of  the  experimental  data.  At 
low  angles  of  attack  the  differences  between  the  GASL  results 
and  the  experimental  data  are  believed  to  be  primarily  caused 
by  viscous  effects.  Since  planes  were  deleted  for  a  >  10  deg, 
comparison  of  leeward  side  pressure  distributions  was  not  pos¬ 
sible  over  much  of  the  body  at  a  =  20  deg  .  Low  pressure  "wells" 
were  found  to  develop  along  the  surface  at  a  >  10  deg  .  The 
experimental  data  did  not  indicate  such  formations;  however, 
further  experimental  data  are  needed  along  the  leeward  side  at 
angles  of  attack. 

5.  The  possible  effects  of  the  low  pressure  "wells"  were  investi¬ 
gated  along  a  constant  entropy  surface  within  the  inviscid  shock 
layer.  The  curvature  of  the  streamlines  and  the  depth  of  the 
low  pressure  wells  were  substantially  less  along  the  constant 
entropy  surface  than  along  the  body  surface.  The  agreement 
between  the  experimental  wall  and  the  computed,  constant 
entropy  surface  pressures  was  better  than  with  the  computed 
inviscid  wall  pressures. 

6.  Force  and  moment  coefficients  were  approximated  reasonably 
well  by  both  the  integrated  GASL  wind -fixed  axis  solution  and 
those  based  on  modified  Newtonian  theory.  The  normal-force  and 
pitching-moment  coefficients,  Cn  and  Cm,  were  found  to  be  in 
good  agreement  with  the  experimental  data  when  computed  from 
the  GASL  results.  In  agreement  with  previous  studies  of 
viscous  effects  on  zero -lift  drag,  the  predicted  inviscid  axial- 
force  coefficient,  ,  was  found  to  be  about  40  percent  below 

the  experimental  data  at  Mn  -  18  over  the  entire  angle -of-attack 
range.  Strong  viscous  effects  were  reflected  in  the  lift-to-drag 
ratio  comparison  where  all  theoretical  results  considerably 
overpredicted  the  lift-to-drag  ratio. 

7.  Normal-force  and  pitching -moment  coefficients,  Cn  and  Cm  , 
were  numerically  differentiated  with  respect  to  angle  of  attack,  a. 
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Comparisons  showed  results  based  on  the  integrated  GASL 
data  were  in  better  agreement  with  the  trends  of  the  experi¬ 
mental  data  than  were  the  results  based  on  modified  Newtonian 
theory. 

The  results  of  the  present  study  indicate  the  GASL  characteristic 
solution  should  be  a  valuable  tool  when  coupled  with  a  three-dimensional 
boundary -layer  theory  for  studying  hypersonic  viscous  effects  at  angles 
of  attack.  Certain  limitations  and  deficiencies  were  noted  which  deserve 
further  study  (such  as  the  comparison  of  pressure  distributions  at  zero 
angle  of  attack). 
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APPENDIX 

FORCE  AND  MOMENT  EQUATIONS 


For  comparison  with  the  experimental  (body-fixed)  pressure  distri¬ 
butions,  it  was  convenient  to  machine  interpolate  the  GASL  wind- fixed 
axis  solution.  The  resulting  surface  pressure  distribution  was  inte¬ 
grated  to  obtain  the  force  and  moment  coefficients  for  comparison  with 
the  experimental  data.  The  transformation,  derivation  of  the  force  and 
moment  equations,  and  the  details  regarding  the  use  of  the  modified 
Newtonian  theory  will  be  given  below. 


TRANSFORMATIONS 

It  can  be  seen  from  Fig.  1  that  the  two  axes  systems  of  interest  are 
related  by  a  rotation  transformation,  which  in  matrix  form  is  as  follows 


lathe  body-axis  system,  the  X  and  V  coordinates  were  replaced  by  the 
polar  coordinates  fin  and  Oq,  so  that 

Xb  =  11b  cos  dB  ,  Yb  =  Rb  sin  Ob  (2) 

where  Xr,  Yb,  and  Rb  are  nondimensionalized  with  respect  to  the  nose 
radius.  Using  the  results  in  Eq.  (2)  in  the  scalar  form  of  Eq.  (1)  the 
result  was  obtained: 


X,v  —  Rg  cos  6b  cos  a  +  Zf)  sin  a  j 

Yft  =*  Yb  =  Rb  sin  0R  '  ( 3) 

Zy,  =  Zb  cos  a  —  Rb  cos  6b  sin  a  ; 

The  geometry  was  such  that 

0B  =  tan-1  (Y'b/Xq),  0.„  =  tan-1  (Yw/Xw) 

Thus, 

(9  iv  =  tan  1  ^Rb  sin  6b >'[ R  b  cos  cos  a  4-  Zp  sina)J  (4) 

so  that  Eqs.  (3)  and  (4)  are  the  transformation  equations  from  body- 
fixed  to  wind-fixed  axes. 

Inasmuch  as  the  reference  axis  for  6W  within  each  change -of -frame 
was  the  shifted  Zw,  there  was  a  correction  on  Xw  in  computing  tan  dw . 
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This  correction  was  to  decrease  Xw  by  the  sum  of  the  AXw’s  for  all 
changes-of-frame  up  to  the  point  considered.  The  result  is 

6vj  =  tan  Rg  sin  0b/(Rb  cos  cos  a  +■  Zg  sin  a  -  2,  AXWl)  (5) 

L  i=  1  J 

where  m  is  the  number  of  changes-of-frame  to  Zw .  For  a  given  Zg 
and  Os,  Eq.  (5)  and  the  last  of  Eq.  (3)  were  used  to  determine  Z *  and  (yw . 
The  pressure  was  computed  by  the  four-point  Lagrangian  formula. 


Y  =  F.WY,  +  F2  (x)  V2  +  Fj(x)  Y3  +  F4(x)Y4 


where 


Ft  (x) 
F,(x) 


.(x  -  x2)  (x  -  x3)  (x  -  x4) 
(x2  -  x2)  {Xj  -  \3)  (x,  -  x4) 

(x  -  x,)  (x  -  x2)  <x  -  x4) 


(xj  x3)  (Xj 


(Xj  -  x4) 


Fj  (x)  = 
F4(x)  - 


(x  -  x2)  (x  -  X3)  U  -  X4) 
(x2  -  X,)  (x2  -  X  j)  (x2  -  x4) 

(x  -  x, )  (x  -  x2)  (x  -  x3) 
(x4  -  x,)  (x4  -  x2)  (x4  -  xs) 


and  the  subscript  numbers  represent  the  known  values  to  be  interpolated 
in  both  Z,v  and  ew .  More  specifically,  for  a  given  Zg  and  0 g  the  corre¬ 
sponding  Zw  and  0W  were  computed.  The  computer  searched  the  GASL 
results  for  four  sets  of  pressures  which  bracketed  both  Zw  and  ,  pref¬ 
erably  twro  on  each  side;  however,  as  the  end  of  the  body  or  a  change- 
of-frame  was  approached,  a  three  and  one  or  an  end  point  interpolation 
was  used.  In  all,  five  interpolations  were  necessary  to  define  the  pres¬ 
sure  at  a  point--four  were  required  to  determine  the  pressure  at  each  of 
the  four  known  6\„’s  in  the  desired  Zw  plane,  and  a  final  interpolation  on  0W 
gave  the  desired  result.  In  the  present  analysis,  these  calculations  were 
performed  at  0.  2  intervals  in  Zg  from  - 1  to  14.644  with  additional  special 
Zb  values  for  comparison  purposes.  In  each  Zg  plane,  pressures  were 
computed  in  10-deg  increments  in  i9g .  The  resultant  pressures  were 
printed  out  and  stored  for  later  use  in  the  force  equations. 


FORCE  COEFFICIENTS 

A  general  expression  for  the  resultant  force  coefficient  attributable 
to  the  surface  pressure  is 

f'F-*  =  S  S  tip  n  d4/jrR|>  (6) 

n  A 

where  rRi/  is  the  reference  base  area,  Cp  the  local  pressure  coeffi¬ 
cient,  and  n  the  local  inward  unit  normal  vector  for  an  element  of 
surface  area,  <IA.  As  can  be  readily  seen,  the  normal  vector  and  the 
elemental  area  are  functions  of  the  local  geometry,  so  that  for  the 
present  case  Eq.  (6)  was  evaluated  for  two  regions --the  spherical  cap 
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and  the  conical  afterbody.  However,  for  ease  of  computing,  three 
regions  were  chosen,  as  shown  in  the  following  sketch. 


In  regions  1  and  2,  the  integration  was  performed  in  the  wind-fixed  axis 
system  because  these  regions  represented  a  sphere  at  zero  angle  of 
attack.  So  that  for  a  given  Zw,  Cp  was  constant,  i.  e. ,  Cp  ?  f  (0*) . 
Immediately,  in  region  1 

2  w, 

(CfZw)  =  -  ( 2/Rb* )  /  Cp  Zw  dZw  (7) 

In  region  1,  the  net  effect  of  the  forces  acting  in  the  Xw  or  Yw  direction 
was  zero. 


The  resulting  equations  for  region  2  are  different  only  because  of 
the  lower  limit  of  integration  formed  by  the  plane  Zb,  ;  that  is,  the  range 
of  integration  over  0W  w^as  +  0W,  to  -0Wj.  For  given  Zw  and  Zb,  ,  one  obtains 


where 


and 


(?B,  =  cos"1  [(Zw  -  Zb,  cos  a)/  -  Rb,  sin  a] 

Rb,  =  -  Zb,1 

=  tan  1  [R  B,  sin  0B,/(Rb,  cos  cos  a  4-  Zb,  sin  a)^ 


Because  of  symmetry,  the  following  equations  result  in  region  2: 


(CFZJ(  =  —  ( 2/  n-  RbS  ) 

(C  F,.}  '  W.Rb*) 


£  Cp  Zw  0W,  dZy 


C  p  \  1  “  Z  yyj  sin 


d  Z  w 


(8) 
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In  region  1,  the  values  of  Cp  were  taken  from  the  blunt  body  (inverse) 
solution,  which  was  used  to  compute  the  initial  data  for  the  GASL  super¬ 
sonic  solution  (region  3).  In  region  2,  the  solution  obtained  from  the 
GASL  program  was  used.  In  region  3,  Cp  =  Cp  ( 6 ,  7,w)  =  Cp  (0r,  Zr>  ,  and 
from  a  geometrical  view  the  equations  are  more  tractable  in  the  body- 
fixed  axis  system.  The  expression  for  the  inward  unit  normal  vector  in 
the  body -fixed  axis  system  at  any  point  on  the  cone  was  substituted  into 
Eq.  (6),  and  the  integration  performed  with  the  result  that 

(C,rz  ^  =  (2  tan  CTc/j7fU2)  f  Rn(Z|})  f  Cp  (^H.  ^d)  JZb 

V  P/3 

(CFxb)  =  -  (2/rrRb’)  J  Rb(Zb)  f  Cp(0u,Zn)  cos  dflodZy 

In  this  region,  the  GASL  results  were  interpolated  for  a  selected  Zb, 
the  integrals  over  evaluated  by  the  trapezoidal  rule,  and  the  results 
were  then  integrated  over  Zr  by  the  trapezoidal  rule.  Because  of  the 
large  number  of  Zq  planes  (150),  this  technique  was  considered  to  be 
sufficiently  accurate.  Because  of  symmetry,  in  all  three  regions  the 
side  forces  integrate  to  zero. 


MOMENT  COEFFICIENT 


The  general  equation  for  the  three  moments  about  a  point  "O"  may 
be  expressed  as 


(Cm)o  -  //  r  x  Cp  n  dA/jrRh*  Db 


(10) 


where  r  is  the  radius  vector  from  point  mO'',  In  regions  1  and  2,  the 
point  "O"  was  chosen  as  the  nose  of  the  sphere  (Zw  =  -1),  so  that 


(Cm)o,. 

(Cm)o,3 


(ID 


The  latter  result  was  obtained  because  the  resultant  moment  arm  of  the 
forces  in  region  2  was  unity.  In  region  3,  the  point  'O"  was  taken  as 
the  intersection  of  Zb  with  the  sphere,  so  that 

B  m  a  x 

(c,n)0  ,  =  -  (2/n-RbSDb)  f  Rb  [Rb  tan  CTC  +  Zy  +  l]* 

ZB.  (12) 

max 

f  Cp  (0B.  Zb)  cos  0r  dZR 


22 


A  E  DC*  T  R-6S-234 


To  add  all  the  contributions  of  the  various  regions,  the  coefficients  in 
regions  1  and  2  were  transformed  to  region  3  through  the  angle  of  attack, 
so  that  the  following  expressions  were  obtained  for  the  axial-force, 
normal-force,  and  pitching-moment  coefficients  and  center  of  pressure 
location  from  the  nose: 


Also 


apd 


Ca 

=  [(Cl?Zv, 

),+ 

(Cfz«L 

|  cos  a  + 

)2  sin  a  +  ( 

A*b), 

) 

Cn 

=  [(CFZw' 

t 

J  v 

),  * 

(Cf  «.),] 

sina-^Cpx 

)2  cos  c  “  ( 

!Cfxb).( 

f 

>  (13) 

Gnlnose 

=  (Cm)o,3 

-[( 

Cfx„X  C0Sa  -  (Cfzw)2 

sine  /  D  h 

\ 

Z ep/Db 

=  Cm/CN 

Cl 

=  Cn  cos  a  -  Ca  sin  a  / 

Cd 

=  C a  cos  O-  +  Cn  sin  a.  * 

(14) 

L/D 

=  c  L  'C  n 

(15) 

These  final  equations  with  the  results  from  the  GASL  program  were 
used  to  compute  the  force  and  moment  coefficients  for  comparison  with 
the  experimental  data. 


NEWTONIAN  THEORY 


The  pressure  coefficient  from  modified  Newtonian  theory  is 

Cp  =  Cp  cos’r;  (  16) 

where  Cp  is  the  pressure  coefficient  at  the  forward  stagnation  point  on 
a  blunt  body.  A  value  of  1. 841  was  used  for  Cp  in  these  computations 
which  corresponds  to  the  normal-shock  pressure  for  M  =  18  ,  The 
angle  v  is  the  acute  angle  between  the  free-stream  velocity  vector  and 
the  local  unit  normal  vector,  i.  e. 


Thus 


cos  r] 


Cp  =  1.841  (V  ■  n  /  |  V  | ) 2 


(17) 


was  used  to  predict  the  surface  pressure  distribution.  The  value  of  Cp 
was  substituted  into  the  force  equations,  and  the  Newtonian  force  coef¬ 
ficients  were  computed.  One  precept  of  the  modified  Newtonian  theory 
is  that  the  pressure  is  zero  in  a  "shadowed''  region.  Thus,  on  inte¬ 
grating  the  pressure  distribution,  the  limit  was  set  at  cos  7  =  o,  or 

Tj  =  77/2  . 
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Experimental  Data 
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Fig.  2  Zero-Lift  Pressure  Distribution  over  9-deg  Half-Angle  Spherically  Blunted  Cones 
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Numerical  Results 
Ideal  Gas  (y  *  1,4),  =  18 

- GASL  Wind-Fixed  Axis  Solution  No.  1 

- GASL  Wind-  Fixed  Axis  Solution  No.  2 

x  >  GASL  Wind-Fixed  Axis  Solution  No.  3 
- GASL  Body-fixed  Axis  Solution 
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Fig.  3  Concluded 
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a  =  10  deg 
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Experimental  Data 


Numerical  Results 
Ideal  Gas  (y  =  1.4),  *  18 

GA5L  Wind-Fixed  Axis  Solution 
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Fig.  5  Pressure  Distributions  along  Most  Windward  and  Leeward  Streamlines  at 
Moc  -  18  and  at  Angles  of  Attack 
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Experimental  Data 


Numerical  Results 
Ideal  Gas  ty  =  1.4),  *  18 

GASL  Wind-Fixed  Axis  Solution 
Modified  Newtonian 
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Fig.  6 


Radial  Pressure  Distributions  over  a  9-deg  Half  Angle 
Spherically  Blunted  Cone  ot  M,-*.-  -  18  and  a  =  5  deg 
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Fig.  7  Radial  Pressure  Distributions  over  a  9-deg  Half-Angle  Spherically 
Blunted  Cone  at  M©o  *  18  and  a  -  10  deg 
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Fig.  8  Radial  Pressure  Distributions  over  o  9-deg  Half-Angle  Spherically  Blunted  Cone  at  M.»  ~ 
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Numerical  Results 
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Fig.  12  Force  and  Moment  Coefficients  for  Spherically  Blunted  Cones  at  M-»  -  18 
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Numerical  Results 
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- Modified  Newtonian 


Experimental  Data 


Sym 

Source 

ReJ\n. 

oc,  deg 

Rn/Rt> 

• 

CAL  48-in.  Shock  Tunnel 

14.7 

14,000 

9 

0.3 

-0.05 

~0.1 

□ 

AEDC-VKF  50-in.  (H) 

18.45 

13,000 

9 

0.3 

-0.05 

-0.1 

AEDC-VKF  100-in.  (FI 

16.5-20 

9,  000-41,  000 

9 

0.3 

-0.05 

-0.1 

▲ 

NASA  Helium  Tunnel 

19 

200, 000 

10 

0.255 

-0.02 

-l.o 

e .  C a  with  a 


Fig-  12  Concluded 
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Experimental  Data 


Numerical  Results 

Sym 

Source 

mgd 

Rero/in. 

Ideal  Gas  <y  -  1. 41,  ,\\x  -  18 

o 

CAL  48-in.  Shock  Tunnel 

14  7 

14.  COO 

Pressure  Distribution  from 

□ 

AEOC-VKF  50-in.  <H) 

18  45 

13,  100 

GASL  Wind-Fixed  Axis  Solution 

Modified  Newtonian 

o 

i 

AEDC-VKF  100-in.  (Ft 
+  10-percenl  Scatter 

16.5-20 

9,  000-41, 000 

Angle  cf  Attack,  a,  deg 


Shaded  areas  are  results  based 
on  experimental  data  with  ±10- 
percent  scatter 


Fig.  13  Variation  of  Normal-Force  and  Pitching-Moment  Coefficient  Slopes  with 

Angle  of  Attack  for  a  9-deg  Holf-Angle  Spherically  Blunted  Cone  at  Moo  '  18 
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TABLE  I 

PRESSURE  DATA  FROM  AEDC-VKF  TUNNELS 


Surface  Distance,  a 

2.  70 

4.  19 

10.  35 

12.  35 

12.35 

12.  35 

12.  35 

14.  28 

15.  28 

10.  35 

Axidl  Distance, 

Zb 

1.  125 

2.  61 

8.  675 

10.  675 

10.  875 

10.  875 

10.  675 

12.  55 

13.54 

8.  675 

Model  Radial  Angle.  e„  .  deg 

0 

180 

0 

157.  5 

112.  5 

67.5 

22.5 

180 

0 

180 

Run  No, 

Tunnel 

limn  el 

Angle  of  Attack 

Roll  Angle 

P/Po 

50-in.  (H) 

136  2* 

0 

0 

0.  0456 

0.  03R0 

0.  0263 

0.  0212 

0.  0223 

0.  0249 

0.  0232 

0. 0227 

0.  0239 

1  3G3 

0 

0 

0.  0459 

0.  0345 

0.  02.38 

0.  0195 

0.  0200 

0.  0221 

0.  0180 

0.  0200 

0.  0183 

13G4* 

0 

0 

0.  0521 

0. 0352 

0.  0249 

0.  0240 

0.  0220 

0.  0240 

0.  0229 

0, 0215 

0.  0232 

1365* 

5 

0 

0. 0769 

0.0212 

0.  0564 

0.0142 

0.  0184 

0.0340 

0.  0513 

0.  0143 

0.  0632 

1366 

5 

0 

0.  0563 

0. 0220 

0.  0459 

0.0113 

0.  01 60 

0.  0312 

0.  0435 

0.  0140 

0.  0542 

1367* 

10 

0 

0.0965 

0.  0213 

0.  1110 

0.  0112 

0.  0147 

0.  0515 

0.  1087 

0. 0110 

0.  1200 

1368 

10 

0 

0.0978 

0,  0184 

0.  1)00 

0.  0102 

0.  0131 

0.  0542 

0.  1016 

0.  0100 

0.  1180 

1369 

5 

0 

— 

0.  0240 

— 

0.  0100 

0.  0143 

0.  0320 

— 

0. 0137 

— 

1370* 

0 

0 

0.0527 

0.  0362 

0.  0236 

0.  0216 

0. 0222 

0.0238 

0, 0221 

0. 0216 

0.  0224 

1372* 

10 

180 

0.0278 

0.  0840 

0. 0113 

0.  1030 

0.  0393 

0.  0140 

0.  0122 

0. 1220 

0.  0110 

1373* 

20 

0 

0.  2100 

0.  0150 

0.  2300 

0.  00G7 

0.  0140 

— 

0.  2200 

0.  0060 

0.  2400 

1374 

20 

0 

0.  2150 

0.  0168 

0.  2540 

0.  0076 

0.  0130 

— 

0.  2310 

0.  0063 

0.  2460 

1379 

10 

30 

0.  0982 

0.  0252 

0. 1070 

0.  0128 

0.  0328 

0.  1060 

0.  1300 

0.0L12 

0.  1160 

1380 

10 

60 

0.  0673 

0.  0269 

0.  0520 

0.  0218 

0.  0710 

0.  1230 

0.  0900 

0.  0129 

0.  0640 

138 1* 

10 

90 

0.  0493 

0.  0407 

0.  0270 

0,  0520 

0,  1080 

0.  1080 

0.  0485 

0.  0300 

0.  0260 

1336* 

20 

180 

... 

0.  2230 

0.  0085 

0.  1950 

— 

--- 

... 

0.  2340 

— 

100-m.(F) 

68b 

0 

0 

0.  0473 

0.  0367 

0. 0234 

0.  0219 

0,  0216 

0.  0235 

690 

0 

0 

0.  0516 

0. 0401 

— 

0.  0233 

0.  0230 

0.  0244 

592 

0 

0 

0. 0494 

0.0389 

0.  0235 

0.  0222 

0.  0220 

0. 0224 

' 

593 

0 

0 

0.  0493 

"  ‘  ~ 

0.0236 

0.  0219 

0.  0215 

0.  0220 

*These  data  presented  in  Ref.  7  also. 
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TABLE  II 

TEST  CONDITIONS  FOR  AEDC-VKF  PRESSURE  DISTRIBUTION  DATA 
(TEST  GAS:  NITROGEN) 


Tunnel 

Run  No. 

Anglo  of  Attack, 
deg 

p0,  psia 

TU,°K 

x  JO*,  p&hi 

rU«  ps|a 

Lr„.  fps 

Pn',  P1*'* 

50-in.  (H) 

1 3  fi  2 

0 

7,  288 

3830 

18.  42 

9,  035 

1. 53 

66,  74 

0.  363 

10,  066 

0.  678 

3  303 

0 

6,  223 

3498 

20.  07 

6,  538 

0.  57 

46.  82 

0.  176 

9,  597 

0.  328 

13G4 

0 

7,  356 

3247 

18.  84 

12,  464 

1.47 

5  3.  37 

0.  366 

9,  206 

0.  682 

1365 

5 

6,  497 

3243 

18.  79 

11,  183 

1.  33 

53.59 

0.  330 

9,  108 

0.  614 

1366 

5 

5,  803 

3460 

21.  76 

5,  701 

0.43 

43.  31 

0.  141 

9,577 

0.  264 

1367 

10 

6,  989 

3897 

18.45 

8,  313 

1.  44 

67.  80 

0.  342 

10,  162 

0.  640 

1368 

10 

6,  937 

4164 

21.57 

4,  526 

0.47 

53.  82 

0.  154 

10,  582 

0.  289 

1369 

5 

6,  609 

4561 

21.  25 

3,  658 

0.  47 

61.  00 

0.  148 

11,  098 

0.  279 

1370 

0 

6,  807 

2862 

19.  41 

14,  274 

1.  23 

44.  05 

0.  324 

8,  614 

0.  602 

1372 

10 

8,  519 

3345 

18.  07 

13,  140 

1.  59 

54.  44 

0.  400 

9,  363 

0.  746 

1373 

20 

7,  263 

3376 

18.  67 

11,  542 

1.51 

56.  65 

0.  367 

9,  400 

0.  685 

1374 

20 

7,  437 

3865 

18.  38 

9,  082 

1.57 

67.  71 

0.  372 

10, 116 

0.  695 

1379 

10 

6,  543 

3473 

18.  99 

0,  320 

1.  19 

56.  47 

0.  300 

9,  547 

0.560 

1380 

10 

6,  519 

2890 

19.  37 

13,  433 

1.  18 

44.  64 

0.  310 

8,  655 

0.  577 

1301 

10 

6,  878 

3315 

18.  80 

11,  U4G 

1.  34 

54.  30 

0.  334 

9,  309 

0.  622 

1336 

20 

7.  523 

3232 

18.  94 

1  2,  680 

1. 46 

52.  58 

0.  366 

9,  185 

0.  682 

100-in.  <F) 

585 

0 

19, 200 

2240 

10.  40 

71,  200 

4.  00 

33.  62 

1.  099 

7,  682 

2.  044 

5  90 

0 

22,  100 

2100 

20.  10 

87,  700 

4.  00 

20.  81 

1.  164 

7,  451 

2.  1G4 

592 

0 

22,  300 

2450 

10.  30 

68,  500 

4.  50 

37.  27 

1.  233 

8,  082 

2.  295 

593 

0 

21,  700 

2480 

19,  20 

66,  600 

4.  50 

38.  22 

1.  235 

8,  112 

2.  300 
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TABLE  III 

MODIFIED  *  PRESSURE  DATA  FROM  CAL  48-IN.  SHOCK  TUNNEL  (REF.  8) 


Surface  Distance,  s 
Axial  Distance,  7 u 
Model  Radial  Auj>lr,  0 

0 

-1.00 

1 . 045 
-0.00 
180 

3.  nlO 

1.  940 
18U 

4.f.25 

.1.025 

0 

4.  G2i 

3.  025 
90 

G  H9I> 

5. 2.7!) 
18!) 

D.  140 

7.  5 00 

0 

9.  140 

7.  500 

90 

11.  390 
9  720 
ISO 

U.  BAH 
11 .  u2.“> 

0 

1  1.  f.  iO 
n.  025 
90 

15.890 

1  1 .  130 
180 

Tunnel 
Kun  iSo. 

Tunnel 

Angle  of  Attack,  deg 

'I  unnol 

Hut  1  Angle,  (leg 

pK  ' 

30 

0 

0 

l  000 

0.  2 3 on 

0.  0437 

0.  0373 

0. 025b 

n,  0345 

0  (1361) 

U  U4J0 

0,  0272 

0  0310 

U  U2lll) 

— 

31 

0 

L 

1. 000 

O.  2olS 

0. O.Wb 

0.  0278 

0. U411 

— 

o.  om 

0  0275 

0.  0 1 H2 

0.  0.(12 

n. 02h  1 

0.  0194 

32 

0 

0 

1 . 000 

o.  rioo 

0.  0204 

0  0245 

— 

0.  0311 

0.  0312 

0.  0234 

0.  024 G 

0  U2U4 

0.  UJsiti 

— 

66 

0 

0 

1.000 

O  24  1(1 

0.  04  H 

0  0293 

0  0392 

0.0127 

0.  03-5 

0.  0257 

0.  0258 

0,  02 9  't 

n.  0*57 

0  0279 

12 

0 

0 

1.000 

0.  24UU 

u.u-isu 

0. 0430 

0.0435 

o  <ur>2 

0. 0415 

0,  0297 

0.  0270 

U.  1)336 

n.  0231  1 

0.  0  307 

13 

5 

0 

0  9920 

(1.  17117 

0.  0771 

0.  0590 

0.  0408 

0.  0’>31 

U. U01J 

0  0292 

O.Olhi 

0  0755 

0  027u 

U.  U1  GO 

14 

10 

0 

0.  0700 

0.  1370 

0.  0240 

0. 0898 

0.  0308 

0.  U18u 

0.  1281 

0.  0290 

0.  0130 

U.  152b 

0.  0297 

0,  0124 

15 

■>o 

0 

0  8830 

O.  079M 

0.0107 

0.  2740 

0  0459 

n. 0154 

0. 3120 

U.  0J98 

0.  0100 

0  3170 

0.U447 

0  OOH.J 

17 

5 

1 80 

0. 0970 

0  280U 

0.  U591 

0.0.725 

0. 0410 

0.  03-it) 

0.  02  In 

0.  0284 

0.  0552 

0. 0205 

0.  02b 5 

0  07  12 

IK 

2U 

180 

0.  88.10 

0.  lb  10 

0  1800 

0.  01B9 

0. 0391 

0.  230j 

0.  U 1  4.1 

0.  (M2  7 

0,  2297 

U.  01U8 

0.  (JOG  1 

(1,  2410 

19 

10 

1BU 

0.0700 

0.  3300 

u  ovso 

0.  0240 

0  0374 

(1.  04, Sh 

0.01  89 

0  02G3 

o.  m» 

0.  01  b7 

0.  0,J7 6 

0.  1343 

21 

10 

99 

0  u?uu 

0.  27hn 

0. OA40 

0.  0409 

_ 

0.  05  30 

U.  U4U) 

... 

0.  CHf.O 

0  0J8G 

U.  U.1K4 

22 

20 

«a 

o.  wno 

0  •’400 

U.  Ob  49 

0.  0519 

--  - 

0.  0802 

0.  0  3  8n 

... 

U  U47b 

0.  05(1  1 

0  0481 

24 

20 

190 

U.  iJK.ICJ 

0,  .'1920 

0. 74?C 

0.  0221 

0  4U10 

U.  01  Jib 

0.  2880 

(J  01  Ob 

0.  2910 

**See  Section  2.  2, 


TABLE  IV 

TEST  CONDITIONS  FOR  CAL  PRESSURE  DISTRIBUTION  DATA  (REF.  8) 

(TEST  GAS:  AIR) 


Tunnel 
Run  No, 

Angle  of 
Attack,  deg 

P's M 

Hu  \  10  \  fri/srri 

M, 

I»m  x  lf|1 .  (■'•u 

"K 

IV  Ps,n 

p,'*.  psu 

30 

U 

2510 

36.  1 

1  1.  65 

11,  583 

2.  84 

13G.  9 

0.  4264 

8401 

0.  7965 

0.  7434 

Hi 

0 

1170 

36.  0 

14.  35 

£>.  fit)  7 

1.  50 

142.0 

0.  2155 

8  381 

0.  4026 

0.  3474 

32 

() 

625 

36.  2 

14.  UO 

3,  083 

0.  02 

150,  1 

0.  1256 

8407 

0.  2346 

0.  2180 

33 

0 

1 177 

35.  y 

14.  37 

5,  667 

1.  40 

141. 1 

0.  2159 

8366 

0.  4033 

0.  3045 

12 

0 

3503 

37.  0 

13.  18 

10.  250 

7.  54 

167.  0 

0.  0587 

8538 

1.  7917 

1. 5261 

13 

5 

3023 

37.  G 

13.  17 

19,  250 

7.  59 

1G7.0 

0.  9641 

855  3 

1 . 8019 

1 . 6090 

14 

10 

3905 

37,  7 

13.  46 

10,  003 

7.  57 

160.  7 

0.  9601 

3569 

1.  7945 

1.  6830 

15 

20 

39  73 

37.  <J 

13.45 

10,  250 

7.  72 

169.  6 

0.  0770 

8584 

1.  8202 

1.  (.350 

17 

ft 

400a 

38.  1 

13.  43 

10,  167 

7.  81 

171. 3 

0.  0863 

8614 

1.  8437 

1. 0510 

18 

20 

4073 

38.  O 

13.  44 

10,  667 

7.  93 

170.  4 

1.  0032 

8599 

1.  8751 

1.  7400 

19 

10 

3737 

37.  5 

13.40 

10,  500 

7.  23 

167.  0 

0.  0101 

8538 

1.  7176 

1.  6170 

21 

10 

4030 

38.  6 

13.  42 

18,  017 

7.  83 

173.  4 

0.  006  7 

8661 

1.  8447 

1.  2950 

22 

20 

3975 

38.  7 

13.41 

18, 500 

7.  74 

174,  3 

0.  0737 

8676 

1.  8205 

1.  3850 

24' 

20 

31190 

38.  1 

13.  45 

18,  583 

7.  52 

17(1.  8 

0.  9518 

8615 

1.  7792 

1.  3730 

TABLE  V 

EXPERIMENTAL  FORCE  COEFFICIENTS  AND  FREE-STREAM  CONDITIONS 


Source* 

T  unnel 
Rim  No. 

AnyLc  oi 
Attack,  dug 

<  A 

Cs 

C.„ 

P|>*  |*s,« 

l'o,°k 

\l„ 

Ri  in. 

x  JO  *  Pb,J 

,  T 

,  (I'-ia 

ji*»r  i 

Pl('* 

CAL 

74 

0 

0,  175 

_ _ 

... 

J192 

36.9* 

11,  75 

13,  655 

3.  69 

138.  1 

0.563 

8495 

1. 051 

100-m.  (F) 

73 

0 

o.  too 

— 

— 

7744 

3650 

10.  1 

0,  ?50 

1.  10 

57 

0.  311 

9800 

0.  580 

uU-m.  (11) 

1223 

5 

0.  195 

0.  HO 

-0.  130 

Gyoo 

3180 

18,  1» 

13,  118 

1. 64 

55 

0.  390 

9115 

0,  727 

CAf. 

68 

5 

0.  190 

0.  114 

-0.  175 

3083 

37.  2* 

14.  73 

14,  25  0 

3,  91 

139.5 

0.894 

8526 

1.  no 

CAT. 

70 

10 

0.  214 

0.  272 

-0.  422 

3800 

37.  G* 

14.  71 

14,  417 

4.  04 

141.4 

0.612 

8:572 

1.143 

100  111.  (F) 

80 

to 

0.  215 

0.  271 

-0.  130 

71 J  7 

3450 

19.  7 

0.  35  0 

1.  00 

52 

0.  273 

9500 

0.510 

100-in,  (F) 

82 

13 

0.  2G 1 

0.  448 

-0.  634 

7y4  3 

3  700 

19.  1 

9.  000 

1.  24 

60 

0,  323 

9900 

0.  b03 

lUU-in.  (K) 

442 

13 

0,  243 

0.  40.4 

-0.  665 

9747 

2566 

16.  !> 

11,  292 

5. 63 

51 

1.  071 

8082 

2.  000 

CAL 

71 

20 

0.  283 

0.  G4.4 

-0.  901 

3592 

3G.  7* 

14.  77 

14,  250 

3,  79 

136.  7 

0.578 

R4fi4 

1.080 

1 00- in.  (F) 

89 

20 

0.  304 

0.  001 

-0.  970 

7474 

3350 

20.0 

10,  450 

1.  OR 

50 

0,  284 

9400 

0.530 

*  It,  w*  nut  given,  ll„  >  10”'  m  IR/ccc'  was.  listed. 


'See  also  Ret.  7  for  AJ2DC-VKF  100-m.  Tunnel  F  anrl  30  in.  Tunnel  H  and  Ref.  8  for  CAL  data. 


TABLE  VI 

COMPUTED  FORCE  AND  MOMENT  COEFFICIENTS  (y  =  1.4) 


Source  of  Surface 
Pressure  Distri¬ 
bution 

Angle  of  Attack, 
deg 

CA 

Cm 

Cm 

L/l.) 

Cn  „/4eg 

Cma/deg 

GASL 

0 

0.  1120 

0 

0 

0 

0. 01620 

-0. 02359 

5 

0.  1182 

0. 0905 

-0. 1337 

0.  635 

0.  02404 

-0.  03731 

10 

0. 1309 

0.  2452 

-0.  3761 

1.  204 

0.  03465 

-0. 05073 

15 

0.  1014 

0.  4284 

-0. 6294 

1.  360 

— 

— 

20 

0.  1976 

0.  6161 

-0, 9075 

1.  290 

— 

- - 

Modified 

Newtonian 

0 

0.  1259 

0 

0 

0 

0. 03042 

-0. 03969 

5 

0.  1308 

0.  1515 

-0.  2159 

0.  972 

0.  02970 

-0.  04128 

10 

0.  1453 

0.  2982 

-0.  4255 

1.  377 

0. 02970 

-0. 04199 

15 

0.  1663 

0.  4520 

-0.  6473 

1. 418 

— 

— 

20 

0.  1901 

0.  6  236 

-0.  8980 

_ 

1.  329 

~  —  ” 

- — 
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